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ABSTRACT 


In this thesis, an optical readout scheme was successfully developed based on the Fourier 
4F optical configuration and integrated with a custom-fabricated microelectromechanical 
system (MEMS)-based, terahertz (THz), detector array. The MEMS THz detector array 
and 4E Eourier optics were able to transduce the THz scene into an optical signal that was 
captured by a commercial charged coupled device (CCD) camera for generating images. 
A quantum cascade laser (QCE) provided the THz illumination for generating images 
while post-image processing performed background subtraction in order to obtain the 
THz scene. The Eourier 4E optical readout system that was implemented was able to 
profile the general shape of the QCE beam pattern and displayed good linearity of 
response of about 23 gray level values per Kelvin. The concept of optical readout from a 
micromechanical sensor array was also validated. 
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I. 


INTRODUCTION 


The terahertz (THz) spectral range, especially within the 0.3 to 10 THz 
frequencies, has continued to generate interest for imaging and spectrometry applications. 
THz wavelengths, which are in the sub-millimeter range, are small enough to be able to 
provide good resolution in imaging and yet also long enough to be able to penetrate thin, 
non-metallic, and non-polarizing dielectric materials like clothing and plastic 
packaging [1]. In addition, numerous explosives and chemical and biological weapons 
have unique THz spectra responses that can be used for detection and identification [2], 
[3]. More importantly, THz radiation is also non-ionizing and potentially safe [4], making 
it attractive for applications seeking human exposure. As a result, there has been 
increasing interest in THz generation and detection. 

To this end, the Sensors Research Lab (SRL) at the Naval Postgraduate School 
has commenced research work in THz sensing and THz imaging. Early research work 
focused on utilizing a quantum cascade laser (QCL) as a THz source and uncooled 
microbolometer cameras as detectors [5], [6]. Subsequently, the research evolved to focus 
on highly absorbing thin films and metamaterial structures. This effort culminated in the 
design and fabrication of novel microelectromechanical systems (MEMS) bi-material 
THz focal plane array (EPA) by Alves et al. [7]-[10]. This EPA consists of an array of 
THz absorbing elements, each attached to bi-material Al/SiOx cantilever legs. These legs 
are anchored to a Si substrate. The absorbing element absorbs the incident THz radiation 
and the heat causes deformation of the bi-material legs. The amount of deformation of the 
legs has been probed by an optical readout system. Similar optical readout techniques 
have also been widely used for IR detection [11]-[13]. 

Based on the custom fabricated EPA, Montagner [14] proposed an optical readout 
design to be integrated with a THz imagining system. The design leverages external THz 
illumination and allows for the displacement of all the sensors in the focal plane array 
(EPA) to be simultaneously measured and to generate an image using a charged couple 
device (CCD) camera. 
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A. OBJECTIVE 

Moving forward, building upon all the work done previously, the aim of this 
thesis is to develop an optical readout scheme utilizing the MEMS based THz sensor to 
perform imaging under THz illumination. 

Specifically, this thesis involves integration of a THz QCL as an illumination 
source, the MEMs sensor, and optical readout system in order to produce an active, real¬ 
time, and standoff THz imaging system operating in transmission mode. 

The proposed setup would specifically consist of a milliwatt scale QCL providing 
active, external illumination onto the MEMS THz EPA. The EPA, which is enclosed in a 
vacuum chamber, would also be illuminated by an external light source and the optical 
system would project the reflected light from the EPA onto a CCD camera. The image 
information from the CCD is then fed onto custom image processing software, designed 
using LabVIEW, to perform image background subtraction for sensor displacement 
visualization and real-time image generation. 

Characterization of the THz imaging system will be done and readout analysis 
will be performed as part of this thesis. 

B. ORGANIZATION OF THESIS 

This thesis is organized as follows. 

• Chapter I provides the background as well as the main objective of this 
thesis. 

• Chapter II provides an overview of different approaches employed in THz 
imaging. 

• Chapter III highlights the main motivations behind the development of the 
MEMS-based THz sensor and the main features of the integrated THz 
imaging system. The various components of the system are also 
introduced. 

• Chapter IV explains the concept of operations of the main components and 
how they contribute to the overall functionality of the integrated system. 
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• Chapter V provides the imaging results and their analysis. 

• Chapter VI presents the eonclusion of this thesis research and suggestions 
for possible future work. 


3 



THIS PAGE INTENTIONALLY LEET BLANK 


4 



II. OVERVIEW OE THZ IMAGING SYSTEMS 


Over the past deeade, THz imaging systems have been developed to provide 
better performanee, simpler operation, and greater portability. This overview provides a 
general survey of THz sourees, THz deteetors and sensors, as well as imaging systems to 
provide an understanding of the important aspeets of THz imaging. 

A. THZ SOURCES 

THz radiation is part of the eontinuous eleetromagnetie wave speetrum. Henee, 
THz ean be emitted by blaekbodies, as deseribed by Planck’s radiation law. THz 
frequencies between 1 to 10 THz are the characteristic peak frequencies of blackbody 
radiation curves of temperatures between 10 to 100 K. Naturally occurring THz sources 
of importance include dust clouds from distant galaxies as well as the cosmic microwave 
background radiation, which is now cooled to below the lower frequency end of the THz 
spectrum. As such, natural THz radiation has important implications in astrophysics [15]. 

For practical imaging applications, THz radiation is usually generated to 
illuminate the object. Some generation techniques include the use of negative differential 
devices, stimulated emission from lasers, acceleration and deceleration of electric 
charges, as well as frequency manipulation by converting other electromagnetic 
frequencies to the THz range. 

The nonlinear current-voltage characteristics of solid-state electronic devices can 
be used to generate THz radiation [16]. THz radiation is generated in this manner by 
Gunn diodes, IMP ATT diodes, and resonant tunneling diodes (RTD). 

THz radiation can also be generated by stimulated emission within a gain medium 

enclosed in a cavity. The generation technique is similar to that developed for lasers 

whereby THz emission takes place when an electron transits from a high-energy to low- 

energy state. The key to successful emission is the maintenance of the population 

inversion in the gain medium by a pump. THz lasers have been generated at low 

temperatures by a Si gain medium doped with P and Bi, and optically pumped by a free 

electron laser [17]. For high-powered applications, optically pumped THz lasers are 
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usually preferred due to their simplicity and high signal-to-noise ratio (SNR). The 
average power for such lasers has been reported to be around 100 mW [18]. 

Accelerating electric charges using magnetic fields can also generate THz 
radiation. Gyrotrons have been shown to produce very powerful coherent and high- 
powered THz beams with peak power between 0.3 to 1.8 kW in 8 ps pulses with a 
frequency range of 0.55 to 1 THz [19]. 

In addition, devices such as frequency multipliers and photomixers can also 
generate THz radiation. Frequency multipliers make use on non-linear devices to convert 
a laser beam into two lower frequency beams. An example of such a non-linear device is 
an optical parametric oscillator such as LiNb 03 , which has been reported to generate THz 
radiation with cryogenic cooling [20]. In the case of photomixers, a pair of optical or 
infrared lasers is irradiated onto a photoconductive semiconductor antenna, usually made 
of GaAs. The laser energies are greater than the band gap of the semiconductor and this 
induces a photocurrent with different frequency. By properly selecting the frequencies of 
the lasers, THz radiation can be generated. The key characteristics of this generation 
technique is that THz radiation can be generated as a continuous wave or as pulses with a 
very long coherence length of up to 8 m at a frequency of approximately 1 THz [21]. 

Quantum cascade lasers (QCLs) are semiconductor-based lasers that can also emit 
THz radiation. Unlike traditional semiconductor lasers whereby photons are emitted 
when electrons from the conduction band recombine with holes in the valance band, 
QCLs exploit electron transitions between sub-energy states within the conduction band. 
This is because quantum confinement splits the conduction band energy levels into sub¬ 
energy levels. This quantum confinement effect is mainly due to the fabrication of 
alternate layers of different semiconductor materials (e.g., AlGaAs and GaAs) with 
different band gap energy. Most QCLs operate based on a three-level energy system. 
Transitions occur sequentially from level 3 to level 1, and the cascading effect is apparent 
when an electron is able to produce many photons as it transits through the levels. The 
limitations of QCL include the need to be cryogenically cooled and a limited frequency- 
tuning range [22]. While QCLs have limited frequency bandwidth, compared to other 
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broad bandwidth systems, QCLs can be tuned to achieve high power for deeper 
transmittance through objects. 

B. THZ DETECTION SCHEMES 

There are two main schemas in THz detection, incoherent or coherent detection. 
In incoherent detection, otherwise known as direct detection, the incident radiation is 
directly converted to heat or used to create charge carriers (i.e., the output from the 
detector is directly proportional to the incident power of the radiation). Examples of 
direct detectors include Schottky diodes, bolometers, thermopiles, pyroelectrics, and 
Golay cells. 

On the other hand, coherent detection involves mixing the incident radiation 
frequency to the frequency from a local oscillator to obtain a lower value, intermediate 
frequency. The intermediate frequency is the difference between the incident frequency 
and the frequency of the local oscillator. Subsequently, the intermediate frequency is 
amplified. Coherent detection not only detects the power of the incident radiation but also 
the phase [23]. 

The incoherent detection scheme is a much simpler approach as compared to the 
coherent detection scheme as it only detects the intensity changes while coherent 
detection captures both the intensity and phase information of the radiation. 

C. CLASSIFICATION AND TYPES OF THZ DETECTORS 

THz detectors can be broadly classified as thermal detectors or photon 
detectors [15]. Thermal detectors generally absorb the incident THz radiation, which 
causes change in some physical property, which can be probed by a readout circuit. 
Bolometers and pyroelectric detectors and Golay cells belong to this category. 

Bolometers are basically made up of materials that change their electrical 
resistance with temperature as they are exposed to radiation. Microbolometers, which can 
be considered as a class of bolometers, are uncooled arrays of bolometers that have 
similar operating principles to those of bolometers. They measure the amount of THz 
radiation by measuring the change in electrical resistance due to temperature change as 
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radiation is absorbed. The most common bolometric materials are made of amorphous 
silicon or vanadium oxides. They generally have a noise equivalent power (NEP) of 
approximately 300 pW in the THz range [24] and 14 pW in the IR range [8]. Hence, 
these microbolometers are generally optimized for IR detection rather than for THz 
detection. A key limitation for microbolometers is that detector FPA has to be integrated 
with a readout circuit, which is used to measure the thermoelectric resistance of the 
detector. This generally makes the fabrication more complex and adds noise to the 
detection. 

Pyroelectric detectors exploit the fact that pyroelectric crystals vary their electric 
polarization as the temperature changes. The pyroelectric crystals are then implemented 
as capacitors with two electrodes mounted perpendicularly to the direction of 
spontaneous polarization. When incident THz radiation is absorbed, there will be a 
temperature change which results in polarization and change in the charge within the 
crystal capacitor. It has been reported that pyroelectric crystals have a responsivity of 1 
kV/W and NEP of approximately 10^ W/Hz'^^ [25]. 

Golay cells consist of a sealed container, containing a gas with low thermal 
conductivity within a flexible window with reflective material. The window is 
illuminated with an external light source and photo detectors track the displacements of 
the flexible window. When THz radiation is absorbed, it heats the gas in the Golay cell 
and causes the gas volume to increase, displacing the window proportionally. The 
reflected light from the window is displaced and the intensity of the incident THz 
radiation can be inferred. For imaging applications, it has been reported that Golay cell 
micro arrays based on a graphene membrane have been developed with a responsivity of 
400nm/K[26]. 

Another class of thermal THz detectors is the MEMS bi-material sensor that 
makes use of meta-material absorbers. This class of detector has been developed by our 
group and will be described in Chapter III. 

The detectors mentioned previously are classified as thermal detectors. The 
second class of THz detectors is photon detectors. Photon detectors produce an electric 
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current when incident THz radiation is absorbed [27]. This oceurs when the ineident THz 
radiation interacts with the eleetrons in the lattice atoms, eleetrons from the dopant atoms 
or free electrons. The incident THz radiation is able to excite these electrons across the 
band gap to induee a photoeurrent [15], which is probed by a readout cireuit. 

An example of a photon deteetor is the Schottky barrier diode, whieh can also be 
used to detect THz radiation. An output voltage is produced in the Sehottky barrier diode 
that is proportional to the power of the input eurrent when it is being indueed by THz 
radiation [28]. It has a non-linear current-voltage relationship, whieh means it can be 
used in a heterodyne reeeiver for eoherent seheme deteetion. In other applieations, it has 
been reported being able to provide wide bandwidth and high sensitivity by eoupling it to 
an antenna and a parabolie mirror to eolleet the THz radiation [29]. 

D. THZ IMAGING 

THz imaging systems ean be broadly elassified as aetive or passive imaging 
systems. Aetive imaging systems illuminate the target with an external THz souree and 
the seattered THz radiation is deteeted. On the other hand, passive imaging relies on 
emitted or seattered THz radiation by a target with no external THz souree illumination. 
Both imaging systems ean operate in either reflection or transmission mode. Generally, 
passive imaging systems rely on room temperature illumination for indoor applieations 
and sky illumination for outdoor applications. 

As a comparison of passive and aetive imaging systems, images from passive 
systems are more influeneed by the object’s thermal emissivity and temperature. Images 
from aetive systems are more influeneed by the angle between the souree of illumination 
and returns received by the detector from the target’s surfaee [30]. Passive imaging 
systems are more suseeptible to thermal noise while aetive imaging systems require high- 
power THz sources to be effeetive. 
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III. IMAGING SYSTEM SETUP 


This chapter introduces the various sub-system components in the proposed 
imaging system setup. The considerations and limitations of each component will also be 
briefly highlighted to derive the rationale for the design of the experimental setup. The 
general features of the system will also be introduced to explain the key concepts of 
operation of the entire system. 

A. MOTIVATIONS BEHIND THE DEVELOPMENT OF THE MEMS THZ 

OPTICAL READOUT SYSTEM 

An alternative to the various detection methods and schemes described in Chapter 
II are MEMS cantilever detectors utilizing optical or capacitive readout schemes. The 
cantilever based MEMS detectors have been successfully developed for infrared (IR) 
detection and imaging [31]. In the optical readout system developed for this thesis, the 
MEMS detector consists of an array of sensing elements (THz absorbers), which absorb 
incoming THz radiation and convert it into heat. The heat is then transmitted by 
conduction to bi-material legs, which undergo bi-metallic deformation when there are 
temperature changes. This bi-metallic deformation is then probed by an optical readout 
by reflecting light off it, as each sensor element has an embedded reflective A1 surface as 
the ground plate. When operated in this manner, each sensing element in the MEMS 
detector acts as a pixel for the optical readout. 

The motivation behind the development of the MEMS THz optical readout 
system is as follows. 

• MEMS-based sensors with optical readout schemes are able to leverage 
commercial-off-the-shelf (COTS) charged coupled devices (CCD), image 
intensifiers (II), or complementary metal oxide semiconductor (CMOS) 
cameras to be used as THz imagers. This can potentially drive down the 
cost of THz imagers and encourage greater adoption. 

• Optical readout schemes allow the MEMS sensor device to be physically 
isolated from the readout system, which is not possible for current 
bolometer based THz imagers. This isolation provides extra design 
freedom and simpler construction, as well as potentially minimizes the 
thermal noise. 
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• Unlike bolometers (which are based on changes in electrical resistance 
with temperature) whereby the responsivity is dependent on the 
temperature sensitive material used in the sensor, the MEMS bi-material 
sensors can be designed with the same materials with different 
responsivity based on the geometrical parameters. 

B. DESIGN OF THE SYSTEM 

The schematics of the overall system and the picture of the actual imaging system 
setup are shown in Figure 1 and Figure 2, respectively. The entire system consists of four 
main parts: 

• The QCF 

• The THz MEMS bi-material FPA 

• Fourier 4F optical readout system 

• Image processing utility software developed using FABview 



Parabolic 
y mirror 1 

i - - >• 

f " 




Parabolic 
mirror 2 


4F optical 
system 


Beam 



Figure 1. Schematic diagram showing the optical readout system. All major 
components are shown here. The transmitted THz beam is absorbed by 
the MEMS detector, which causes thermo-mechanical deflection of the 
bi-material sensor. This deflection is probed by the Fourier 4F optical 
system and fed to the CCD camera for imaging. 
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Figure 2. Actual imaging system setup, which was implemented in the 
laboratory. The Fourier 4F optical readout assembly is offset to 
compensate for the intrinsic tilting of the absorbers due to residual 
stress during fabrication. 


C. GENERAL FEATURES OF THE SYSTEM 

The following general system features provide real-time and high-frame rate 
imaging capabilities. 

1. Optical Readout from Sensor 

Optical readout from the sensor enables it to maintain high sensitivity with very 
low crosstalk. Optical readout schemes are also relatively simple to design as compared 
to complex, integrated microelectronic readout schemes. In this setup, the images 
obtained from the readout are fed to image-processing software for downstream 
processing and image enhancements. 


13 












2. Multiple Pixel Probing 

The MEMS sensor is physieally isolated from the optical readout system and the 
displacement of each pixel is optically tracked. The optical readout system is able to 
process the displacement of each pixel simultaneously. Hence, multiple pixels could be 
probed at the same time in order to generate the THz scene. 

3. Differential Imaging 

Due to the use of bi-materials, the sensors in the MEMS detector have residual 
stress and this causes the legs and the absorber of the sensor to have an intrinsic 
deflection. As the sensors absorb THz radiation, the bi-material legs increase in 
temperature causing further deflections. Differential imaging was used to eliminate the 
intrinsic bending pattern and detect changes from the initial resting position of the sensor. 

Eor the setup, the differential imaging involves capturing the background 
information initially and storing it in memory. This background image is then subtracted 
from subsequent images captured in order to produce a difference between the two. This 
difference reveals the displacement information of each sensor on the THz EPA. In turn, 
the displacement information directly relates to the amount of THz radiation, which is 
absorbed by the sensors at each location. This concept of operation is further elaborated 
in the later part of this chapter. 

D. COMPONENTS OF THE SYSTEM 

The individual components used in system are highlighted in the following 
sections. 

1. The Sensor 

The sensor used in the system is the heart of the THz imaging system. The MEMS 
bio material EPA used in this work was designed and fabricated by Alves et al. [32]. It is a 
thermo-mechanical sensor, which converts THz energy absorbed into mechanical 
displacement. The performance parameters of the sensor are strongly influenced by the 
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materials used in the detector and the geometry of the design. These performance 
parameters are shown in Table 1. 


Table 1. Measured performance parameters of MEMS detector. From [32]. 


Absorbance, t| 

Time 

Responsivity, 

Thermo-mechanical 

Noise Equivalent 


constant, 

dO/dP 

sensitivity, dT/dP 

Power, NEP 


T (S) 

(X 10® Deg W') 

(Deg K ') 

(X lO-^W) 

0.95 

0.8 

1.2 

0.2 

8.6 


The MEMS detector used consists of an array of THz absorbers. Each absorber 
element also has 10x10 square meta-material sensing elements, which were optimized 
to absorb at 3.8 THz of the QCE. In addition, each absorber is connected to the Si 
substrate by rectangular bi-material legs to a supporting anchor made of SiOx. The bi¬ 
material legs and the anchor conduct heat away from the absorber to the Si substrate. The 
Si substrate acts as a heat sink for this purpose. Incident THz radiation is converted into 
heat by the absorbers. This heat is then transmitted through the bi-material legs. The 
transmitted heat causes the temperature of the legs to increase and undergo bi-metallic 
deformation. This, in turn, alters the angular deflection of the absorber, which is then 
probed by an optical readout by reflecting light by the A1 ground plane. The structure of a 
sensor in the FPA is illustrated in Figure 3. 
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Figure 3. MEMS bi-material THz sensor pixel used for the system (all 
dimensions are in pm), (a) Top view of each sensor, (b) Vertical cross 
section of the absorber and bi-material legs. Square meta-material 
structure on top of the absorber provides strong THz absorbance while 
the Al ground plate provides the reflective surface for optical readout. 

From [32]. 


The performance parameters highlighted in Table 1 are independent from the 
optical readout except the NEP. It establishes the intrinsic limits of the sensor, which are 
dependent on the structure of the sensor and the materials used. The intrinsic noise of the 
sensor includes temperature fluctuation noise, background fluctuation noise, and thermo¬ 
mechanical noise [33]. Temperature fluctuation noise is a result of conductive heat 
exchange of between the sensor and the surroundings due to the temperature fluctuations 
between them. As a result, the sensor tip constantly fluctuates. Thermo-mechanical noise 
is an inherent property of the sensor caused by the exchange of mechanical energy 
between the sensor and the external environment. In this case, the mechanical structure 
acts as a resonator, which can store and dissipate mechanical energy. The oscillations of 

the cantilever structure are directly proportional to the thermal energy of the cantilever. 

16 




























Background fluctuation noise is caused by the temporal fluctuations in the radiative heat 
absorbed by the sensor. The primary heat exehange mechanism in this case is radiation 
between the sensor and the background. Temperature fluctuation noise, thermo- 
mechanieal noise, and background fluctuation noise contribute to the overall fluetuations 
in the sensor deflection. This differs from the readout noise, which is mainly due to 
illumination source fluctuation and mechanical vibrations in the readout system. 

2. Sensor Vacuum Chamber 

Temperature fluctuation noise, background fluctuation noise, and thermo¬ 
mechanical noise cannot be prevented. However, it is still possible to mitigate convective 
heat exchange between the sensor and the ambient surroundings. To mitigate the effeets 
of convective heat exchange, the sensor is placed in a vacuum chamber at a pressure of 
approximately 4 x 10-^ mbar in order to minimize heat exchange between the sensor and 
the surrounding ambient environment. The vacuum chamber was custom fabricated at the 
Naval Postgraduate School and is shown in Figure 4. 



To vacuum 
pump 


Tsurupica 

window 

(behind) 


Quartz window 


(front) 



(a) 


(b) 


Figure 4. Vacuum chamber used to house sensors to minimize convective heat 
exchange, (a) Closeup front view of vacuum chamber with mounted 
MEMS THz sensor array with the front and back transmission 
windows, (b) Vacuum chamber mounted on precision rotating stage for 

alignment purposes. 
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The vacuum chamber has a front tsurupica window and back quartz window. The 
front tsurupica window has approximately 80 percent transmittance at 3.8 THz [34] to 
allow most of the QCL THz radiation to reach the MEMS detector. The quartz window 
has an anti-reflective coating to allow maximum transmission of visible light from the 
readout side of the detector. Finally, the entire vacuum chamber is mounted on a 3-axis 
micrometer driven rotating stage to facilitate alignment as shown in Figure 4(b). 

3. THz Source - Quantum Cascade Laser (QCL) 

A THz-QCL is used as the illumination source since it provides the most suitable 
means to generate THz radiation. QCF-based THz sources are simpler and more compact 
than other sources (see Figure 5) and can provide power at a milliwatt range. 



Frequency (THz) 

Figure 5. Plot of output power against THz frequency for the various methods 
to produce THz radiation. It can be seen that from the required output 
and frequency, THz-QCF provides the most suitable means. From [35]. 

A QCF of approximately 3.8 THz is used as an illuminating source during the 
measurements. The continuous electrical pumping of the QCL results in rapid joule 
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heating, which causes degradation of output power. Therefore, the QCL is operated in 
pulsed mode and cooled down to approximately 10 to 15 K for effective operation. 
Figure 6 shows the QCL and the cryostat assembly. 



Figure 6. QCL and cryostat assembly. The QCL needs to be housed in the 
cryostat chamber to be cooled for effective power output. The QCL 
head is visible through the tsurupica window of the cryostat assembly. 
The tsurupica window is used to allow maximum transmission of the 
generated THz radiation from the QCL. 


The QCL is driven by an AVTECH pulsed laser diode driver, AVO-6F1Z-B. It is 
also gated with an Agilent Technologies 33522A function/arbitrary waveform generator. 
For normal operations, the QCL is operated at a pulse repetition frequency (PRF) 
between 0.5 to 10 kHz, at a constant pulse width of 5 ps and gated at 500 mHz to prevent 
the laser from overheating so that the laser performance does not degrade. However, for 
short durations of approximately 15 to 20 seconds, the gating can be switched off for 
generating greater power. The QCL performance parameters used are summarized in 
Table 2. 
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Table 2. QCL parameters used for THz imaging. The measured peak 
frequency and spectral width were obtained from [36]. 


Peak 

frequency 

(THz) 

Full width 
half 

maximum 
spectral 
width of 

QCL (GHz) 

Pulse 

repetition 

frequency 

(kHz) 

Operating 
laser Pulse 
width (ps) 

Output 
power (pW) 

QCL gating 

frequency 

(mHz) 

3.78 

30 

0.5-10 

5 

~1 

500 


4. THz Beam Focusing Optics 

The QCL source generates THz radiation as a highly divergent beam in a 
concentric ring pattern as shown in Figure 7 (from [36]). 



Figure 7. Imaging of QCL THz laser beam using a microbolometer camera. 
The beam was imaged using two converging lenses by Ng [36]. It can 
be observed that the THz beam pattern consists of concentric rings. 
These rings are caused by far field interference due to multiple facets of 

the QCL. 
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To focus the THz beam energy onto the sensor, it is neeessary to use foeusing 
opties. Off-axis parabolie (OAP) mirrors were seleeted beeause they are able to refleet 
light from a point souree and convert it to a eollimated beam off axis. The eollimated 
region is to provide an area for the target to be imaged. The OAP mirrors used are gold 
plated for maximum reflectanee. 

For the setup, two OAP metal mirrors (fl and f2) are used to foeus and steer the 
beam. The fl (foeal length=50.8 mm) OAP is used to eollimate the emerging THz beam 
while the f2 (foeal length=101.8 mm) OAP is used to foeus the beam onto the sensor. 
The OAP mirrors are configured to be as elose as possible with each other to minimize 
the optical path so as to reduce beam losses due to absorption. The use of OAP mirrors is 
preferred over lenses to focus the THz beam beeause the OAP mirrors absorb less THz 
radiation as compared to tsurupica lenses. The downside is that OAP mirrors are much 
harder to align. 

5. 4F Optical System 

Figure 8 shows the optieal readout system used for the THz imaging. It is 
eomposed of a eollimated light souree, a beam splitter, a pair of lenses (foeal length=10 
em), and an aperture, whieh acts as a spatial filter. 



Figure 8. Optieal readout system. The system is mounted on an optical rail for 

alignment purposes. 
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The 4F optical system is a “perfect imaging system,” whereby all rays emerging 
from a single point in an object space would re-converge and form another point in the 
image space after passing through the system [37]. This allows the optical scene from the 
sensor (input plane) to be imaged onto the CCD camera sensor (output plane). Moreover, 
the Fourier plane between the lenses allows spatial filtering to be performed. Spatial 
filtering allows the system to have increased sensitivity for the optical readout. 

For our application, an aperture is used in the Fourier plane to serve two purposes. 

• The aperture would be able to block deflected rays from the sensor and 
thus provide spatial filtering for enhanced sensitivity. This concept of 
operation will be further elaborated in Chapter IV. 

• The aperture acts as a low pass filter and applies a Gaussian blur to the 
image to reduce image noise [38]. This is a simple approach to reduce 
image noise without having to implement digital processing downstream 
in the system. In addition, it provides better control for the user to 
determine the amount of detail to be revealed during the imaging. 

The main limitation of utilizing a 4F optical system was that it was relatively 
difficult to align it with the MEMS THz sensor. 

6. Collimated Optical Light Source 

The light source used for readout is a high-power light-emitting diode (LED), 
CREE 941-XPEGRNL lOOAOl, biased at a voltage of approximately 3 V and current 350 
mA. This allows the LED to emit green light at a theoretical maximum of 87.4 Im [39]. 
The LED is placed behind a small aperture to act as a point source. The light emitted is 
collimated by a converging lens before it reaches the beamsplitter. The collimated LED 
assembly is shown in Eigure 9. 
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Figure 9. Top view of the self-aligned eollimated LED assembly. The LED 
point souree was enelosed in a lens tube produee a eollimated output. 
The beam splitter refleets the eollimated light onto the MEMS deteetor 
and the refleeted light from the MEMS deteetor is then transmitted 
through the beamsplitter to the Eourier 4E optieal system. 


7. CCD Camera 

The model of the CCD camera used is a Easier A631f The purpose of the CCD 
camera is to convert the image into a digital format for post processing. The CCD camera 
has 1392 x 1040 pixels. Each pixel measures 4.65 x 4.65 pm, which is of sufficient 
resolution to image each pixel from the MEMS THz sensor EPA. The output of the 
camera has 8 bits per pixel, which allows a maximum of 255 gray values. These gray 
values allow the indirect measurement of the relative intensities of the reflected light on 
the CCD sensor. In addition, the camera has a maximum frame rate of 18.7 fps, which is 
fast enough for our imaging application. To complement the maximum sensitivity of the 
CCD camera, the green light source from the EED is used to illuminate the sensor. 
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Differential Imaging Software 


The digital signal from the CCD camera is fed to custom-designed imaging 
software developed by Montagner [14], The imaging software is programmed using 
LabVIEW, developed by National Instruments. The high-level flowchart describing 
differential imaging performed in LabVIEW is shown in Figure 10. 



Figure 10. Steps used in differential imaging using LabVIEW. Step I requires 
the user to determine the gain and shutter settings of the camera. These 
settings affect the image quality. Step 2 obtains the background image 
without illumination and step 3 obtains the image with illumination. 
Step 4 does the most important function to form the THz image. Step 5 
involves detecting local discontinuities in the pixel intensities by 
allowing the user to determine the range of pixel gray values that would 
be replaced by a specific one. From [14]. 
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The imaging software works by capturing an initial frame as a background image. 
Subsequent images acquired by the camera are processed by subtracting from the 
background image. This is an indirect means to visualize the displacement of the sensors 
in the FPA, which is a result of the THz scene on the FPA. For the subtracted image 
display, there are two options available: 

• Simple subtraction whereby the gray scale value recorded by each pixel in 
the CCD sensor is subtracted from the gray scale value from the same 
pixel in the background image acquisition. The result can either be a 
positive or negative value. 

• Absolute subtraction whereby the gray scale value is the absolute value 
obtained from simple subtraction. 

The developed software could also perform image edge detection in the final step. 
The image edge detection could identify edges by associating big changes in pixel 
intensity values as edges. For the setup, absolute subtraction is preferred as it is the 
simplest means to measure the absolute displacement of the sensor, which corresponds to 
the relative intensity of THz radiation received by the sensors. This is a suitable means to 
visualize the amount of THz radiation incident on the sensor, as the setup would mainly 
operate in transmission mode. 

In addition, the imaging software allows a specific area on the obtained image to 
be surveyed. This allows the statistical measurement of pixel intensities of the obtained 
image for characterization purposes. The graphical user interface of the imaging software 
is shown in Figure 11, and a screenshot showing a specific area survey for statistical 
measurement is shown in Figure 12. 
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Figure 11. Graphical user interface of imaging software developed using 
Lab VIEW. The results from the absolute subtraction due to sensor 
deflection are shown in the last screen on the second row. Various 
functions are available to allow the user to adjust the image quality 

obtained. 


Absolute Subtraction 



11392x1040 IX 8-bit image 0 (1106,7^ 


Figure 12. Screenshot showing a specific area (selected in the green box) for 
statistical measurement. The mean intensity and standard deviation of 
the gray values within the green box will be calculated by the program 

and provided to the user. 
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IV. IMPLEMENTATION AND CONCEPT OE OPERATION OF 

THE IMAGING SYSTEM 


This chapter details the specific measures used to integrate the various 
components of the imaging system and the concept of operation for the imaging system. 

A. FPA AND OPTICAL SYSTEM ALIGNMENT 

The imaging system is implemented with an angular offset between the vacuum 
chamber housing the sensor and the 4F optical readout system due to intrinsic bending 
(due to residual stress) of the bi-material pixels. The intrinsic bending is measured and 
found to be approximately 8° with respect to the substrate plane [32]. This is shown in 
Figure 13. 



Figure 13. 3D profile plot of sensor showing the intrinsic bending of the sensor. 

From [9]. 

The sensor FPA is positioned perpendicular to the incident THz radiation and the 
optical readout system is offset at an angle. This differs from another approach adopted 
by Grbovic [33] to compensate for the intrinsic deformation of an IR sensor during 
fabrication. In that approach, the angular position of the FPA was adjusted instead. This 
difference is illustrated in Figure 14. Grbovic’s approach is illustrated in Figure 14(a) 
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while the approaeh adopted for this thesis is illustrated in Figure 14(b). The offset in the 
aetual setup ean be seen in Figure 15. 



(a) 


(b) 


Figure 14. (a) Illustration of FPA angular adjustment adopted by Grbovic to 

eompensate for intrin s ic stress. From [33]. (b) Illustration of FPA 
angular adjustment adopted. 
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Figure 15. Top view of aetual setup of optieal system implemented with 

offsetting opties with the FPA. This is done to ensure maximum THz 
radiation ineident on the sensor and inerease sensitivity of the system. 


The rationale for the angular offset of the optieal system instead of the sensor 
FPA is as follows: 

• To ensure maximum amount of THz radiation ineident on the sensor. 

• To maximize refleetion from the ground planes of the pixels and henee 
small intensity ehanges eould be deteeted easily. 

As shown in Figure 14, the best image quality is obtained when the offset angle is 
10.7°. This value is elose to the measured bending of the sensor of about 8°. 
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B. CONCEPT OF OPERATION OF 4F OPTICAL SYSTEM 


The concept of operation of the 4F optical system is illustrated in Figure 16. 
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Figure 16. Concept of operation of 4F optical system. The first lens performs 
Fourier transform of the object and generates the spatial frequencies at 
the Fourier plane. The second lens performs inverse Fourier transform 
to reconstruct the original image. Spatial filtering is done at the Fourier 
plane. Due to deflections of the sensor, reflected rays are focused onto 
different positions on the Fourier plane. This allows spatial filtering to 

be performed. 


The reflected light from each sensor on the FPA is decomposed into its respective 
spatial frequencies after passing through the first lens. Essentially, the first lens performs 
Fourier transform to produce the spatial frequencies of the object to be imaged. This is 
expressed in Equation (1) 

f(«, v) = (1) 
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where F(u,v) is the image function, u and v are the two dimensional coordinates in the 
output plane of the first lens, 2 Lndf(x,y) is the two dimensional coordinates in input plane 
of the first lens. 

Subsequently, the spatial frequencies are focused at the Fourier plane by the first 
lens and the second lens reproduces the image at the output plane (FPA of the CCD). The 
second lens basically performs another Fourier transform operation, which essentially is 
the inverse Fourier transform of the spectrum at the Fourier plane. This produces the 
original image at the output plane. 

By regarding each pixel on the detector as a point source, each pixel when imaged 
through the optical system is an impulse response. The output signal is then a summation 
of all the impulse responses of the pixels. Mathematically, it is a convolution of the input 
signal with an impulse response as in Equation (2): 

output{x,y) - input{x,y) 0 h{x,y) (2) 

where output(x,y) is image function at the output plane of the 4F optical system, 
input(x,y) is the object function at the input place, and h(x,y) is the impulse response of 
the optical system. In the proposed design, a 2-axis circular aperture filter with adjustable 
aperture size is used to allow selective filtering of frequencies. The impulse response of 
the circular filter is basically an Airy disc and is described in Equation (3) from [40]: 
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( 3 ) 


2 2 

where r = x + y , Ji is Bessel function of the first order, D is the diameter of the aperture 
filter,/is the focal length of the lens used, and X is the wavelength of the light used. 

Based on the convolution theorem, the Eourier transform of a convolution of two 
functions is just the multiplication of the Eourier transform of the functions in the spatial 
frequency domain. In Eourier space, 

OUTPUT {u, v) = INPUT {u, v)H(u, v) (4) 
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where OUTPUT(u,v) is the Fourier transform of the output image, INPUT(u,v) is the 
Fourier transform of the image speetrum at the input plane, H(u,v) is the Fourier 
transform of the impulse response, otherwise known as the eoherent transfer funetion. For 
the Fourier 4F readout system in this thesis, H(u,v) can be described as the filtering 
function of the optical system. In this case, the filtering function changes the value of 
specific spatial frequencies to be zero. 

With the mathematical descriptions in place, it is appropriate to provide the 
concept of how the Fourier 4F optical system performs spatial filtering. From Figure 18, 
the concept of the detector image being converted into its respective spatial frequencies at 
the Fourier plane is shown. 




Figure 17. Illustration of an image from the detector decomposed into its 

respective spatial frequencies, (a) Original image of light reflecting of 
the pixels from the sensor before passing through the first lens for 
Fourier transform operations, (b) Fourier transform of the image 
obtained after passing through the first lens. Highlighted are the relative 
locations of the spatial frequencies. The Fourier transformation of the 
image was obtained using MATLAB. 
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At the Fourier plane, the low spatial frequeneies are mainly due to the eentral part 
of the sensor, whereby most light is reflected off the absorber ground plane while high 
spatial frequencies are due to the edges of the absorber. The effect of high-pass and low- 
pass filtering can be seen from the FPA images in Figure 18. It is apparent that the 
presence of higher frequencies reveals more details in the image, especially at the edges 
of each sensor. In addition, a larger aperture also allows more light to reach the CCD. 
Hence, Figure 18(b) is a brighter image compared to Figure 18(a), even though the 
amount of illumination remains the same. 



Figure 18. Reconstructed images of the same original obtained after passing 
through the Fourier 4F optical readout system with different aperture 
sizes, (a) Low-pass image obtained with a small aperture opening, (b) 
Broadband-pass (high + low pass) image obtained with a large aperture 

opening. 


More importantly, for sensitivity in change detection, the aperture diameter must 
be small enough. Changes in the position of the sensor would cause the light rays from 
the sensor to be focused at different positions at the Fourier plane. If it hits opaque 
positions on the filter, it will not be reproduced at the output. Hence, an aperture is used 
to partially block rays due to deformations in the sensor, which gives the gray scale. 
Small apertures increase the sensitivity of the optical system at the expense of reducing 
the amount of light for image reproduction at the CCD camera. 
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To calculate the theoretieal sensitivity of the system, starting from first prineiples 
of opties and regarding eaeh pixel as a “pseudo slit,” the first minimum is given by 

e^- (5) 

a 

where 6 is angle between the normal line of the sereen and ray from the pixel to the 
image on the Fourier plane, X is the wavelength of the light from the eollimated LED 
assembly, and a is the width of the pixel. 

When eonsidering the first lens in the optieal system, the angular width of the first 
minimum on the Fourier plane then beeomes 


a 

where/is the foeal length of the first lens. 


( 6 ) 


Using Rayleigh’s eriterion, to be able to resolve the reflected light displacement 
from the pixel properly, the zero order maximum must shift to the first minimum of the 
original image. Using Equation (6), with the estimated parameters of i = 500 nm, f = 10 
em, and a = 200 pm, the optical system is estimated to be sensitive enough able to deteet 
the pixel’s angle of defieetion of approximately 1.43 x iQ- deg, whieh eorresponds to a 
temperature ehange of 2.86 x 10- K on the surfaee of the sensor. 

In summary, the 4F optieal system is a preferred optieal arrangement for the 
following reasons. 

• Spatial filtering ean be performed at the Fourier plane. This makes 
downstream imaging proeessing simpler. 

• There is no need for foeusing opties for the CCD eamera as the output rays 
are collimated. 

• The sensitivity of the system ean be easily adjusted by varying the 
diameter of the aperture filter. 
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C. OPTICS ALIGNMENT PROCEDURE 

To ensure the proper funetioning of the system, various optieal eomponents are 
required to be aligned and eonfigured in the following sequenee. 

1. ALIGNMENT OF QCL BEAM WITH FOCUSING MIRRORS 

As the THz imaging system is operated in transmission mode, it is neeessary to 
produee a broad, eollimated THz beam to illuminate the target and produce the THz 
scene on the MEMS THz sensor. This requires precise alignment of the off axis parabolic 
(OAP) mirrors. 

The alignment of the QCL beam with the focusing OAP mirrors is carried out as 
follows: 

• The fl (focal length=50.8 mm) OAP mirror is placed at the exit window 
of the QCL. A semiconductor laser affixed to a beam broadener provides a 
highly collimated light source. The collimated light source is directed on 
the f/1 OAP mirror and adjusted in such a way that the focal point of the 
OAP is formed on the QCL. In addition, to ensure the maximum amount 
of THz radiation is captured by the fl OAP mirror, the collimated light 
needs to illuminate near the center of the OAP. 

• Next, the f2 (focal length=101.8 mm) OAP mirror is placed approximately 
10 cm away from the fl mirror. The collimated light source is then 
directed on the f2 OAP mirror and adjusted such that the focal point of the 
f2 mirror is near the center of the fl mirror. 

• A microbolometer with Ge lens removed, exposing the microbolometer’s 
sensor LPA, is used to image and confirm the QCL beam pattern produced 
from the alignment. 

2. ALIGNMENT OF 4F OPTICAL SYSTEM 

The alignment of the 4L optical system is carried out as follows: 

• The 4L optical system is initially placed perpendicular to the MEMS THz 
sensor. The system is then focused on the substrate plane of the sensor. 

• Subsequently, the entire 4L optical system is rotated until the pixels from 
the sensor appear on the image. Line adjustments are made and the 
aperture filter is closed as small as possible. 
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• The QCL is then turned on and additional adjustments are made to 
aperture size and position so as to obtain the greatest sensitivity. 

D. CCD CAMERA SETTINGS 

The camera setup parameters are configured through the camera setup options in 
the imaging software. The shutter value on the CCD camera is set to a maximum value of 
4095. The exposure time is calculated based on the following formula [41]; 

Exposure Time = (shutter value setting) x 20 ps (7) 

Using Equation (7), this translates to an exposure time of 81.9 ms or 12.18 fps. 
Eor the setup, the maximum camera exposure time is used to collect more light to form a 
brighter image for better differential imaging. Moreover, at this setting, the speed of 
operation of the CCD is still much faster than that compared to the MEMS THz sensor 
(time constant, x = 0.3 s) [7]. Thus, the CCD camera operating at the exposure limit is 
still fast enough for real-time imaging using the MEMS THz sensor. 

The camera’s sensor has 10-bit dynamic range while the camera output has an 8- 
bit output range. As the diameter of the aperture used is small (up to approximately 0.5 
mm), the amount of light reaching the camera CCD is also correspondingly very low. In 
order to increase the detectability of light, it is necessary to map a smaller output of the 
camera’s sensor output range to the full dynamic range of the camera output. Thus, a 255 
gain value is applied. This corresponds to a 12 dB gain (a 4x amplification) and the 
output response is shown in Eigure 19. 
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8 Bit Camera Output 



Increasing Light 


Figure 19. Output due to gain value of 255 being applied. From the graph, it 
ean be observed that the eamera output is fully saturated when the 
sensor’s output is only 25 percent. This produces an effective 12 dB 
gain or 4x magnification. From [42]. 
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V. RESULTS AND ANALYSIS 


This chapter outlines the test procedures and presents the results obtained from 
the integrated imaging system. An analysis of the results will also be given to draw 
important conelusions regarding the eharacteristics and performanee of the integrated 
imaging system. While the intrinsic properties of the MEMS THz sensor have been 
presented in Chapter III, the focus of the eharacterization is to assess the performance of 
the integrated optical readout system. 

A. IMAGING RESULTS 

Using the integrated optical readout system, imaging is performed while 
subtracting the background to eliminate the fixed pattern noise. Figure 20(a) shows the 
image of a QCL beam captured from a commercial microbolometer camera and Figure 
20(b) shows the image of the similar QCL beam captured by our integrated optical 
readout system, whieh was operated at a PRF of 5 kHz and gated at 500 mHz. 




(a) (b) 

Figure 20. Imaging of QCL beam. The diffraction pattern is due to far field 
interference caused by the facets of the QCL laser. The QCL beam is 
adjusted to give a ‘mushroom’ shape in order to be imaged on the 
MEMS THz detector, (a) Image obtained using a eommereial IR 
microbolometer camera, (b) Image obtained using MEMS THz deteetor 
with background subtraction. It can be observed that, due to the larger 
pitch size of the MEMS THz FPA as compared to the microbolometer 
FPA, it eannot resolve the finer details of the QCL beam. 
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Due to the larger spacing between the pixels in the MEMS THz sensor, the 
resolution is much lower than that compared to the QCL beam pattern imaged by the 
microbolometer camera. An attempt is also made to image objects under THz 
illumination. As THz radiation cannot penetrate through metal, a thin metal rod is placed 
between the OAP mirrors to act as a target (see Figure 21). The metal rod (diameter ~ 2 
mm) is then moved about to be imaged on the detector and incite a response. The gating 
of the QCL is switched off for the QCL to operate in continuous mode to see the 
movement of the metal rod more easily. Figure 21 shows the results of the imaging for 
the moving metal rod. 



Figure 21. A thin metal rod moving about and imaged under THz illumination 
using the integrated optical readout system. The position of the metal 
rod is highlighted in yellow, as the resolution of the system is not able 
to provide a clearer outline of the metal rod. As the metal rod moves 
about the illumination scene, the pixels in the detector start to flicker. 
The ‘flickering’ effect is due to the metal rod blocking THz radiation. 

This effect is more apparent in the video recording of the frames 

captured above. 


From Figure 21, it can be observed that the effects of the large pitch of the sensor 
limit the resolution of the integrated optical readout system. The system is unable to 
resolve the outline of the metal rod. However, as the metal rod moves through the 
illumination scene, it blocks THz radiation and causes the sensors to deflect. This makes 
the THz seem to ‘flicker.’ Hence, while the optical readout system suffers from poor 
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resolution, it is sensitive enough to detect small changes in THz radiation due to a small 
moving metallic object. 

B. SYSTEM THERMAL RESPONSE CHARACTERIZATION 

The thermo-mechanical sensitivity (dT/d0=O.2 Deg K'*) of the MEMS THz 
sensor was presented in earlier chapters. It is an indication of the thermal sensitivity of 
the pixel itself and does not fully describe the thermal sensitivity of the integrated readout 
system. It was expected that the integrated thermal responsivity of the entire system, if 
integrated and implemented well, would approximate to the value obtained for the 
thermo-mechanical sensitivity of the MEMS THz sensor. However, instead of monitoring 
the angle of deformation of the sensor as the temperature varies, the change in the grey 
pixel value is used instead. The thermal response characterization aims to measure the 
change in grey level values per degree change in temperature of the sensor. The value 
obtained is the combined thermal sensitivity of the MEMS THz sensor and the 4E optical 
readout system. To obtain the thermal response of the system, the setup is configured as 
shown in Eigure 22. 

4f Optical 
system 


Beam 



Eigure 22. Imaging setup for thermal response characterization. 4E optical 

system is used to detect changes in the sensor deflection as temperature 

changes. 
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Figure 23. MEMS detector mounted on copper heating element to monitor 
thermal response of sensor in 4F optical system. 

The temperature of the sensor can be varied by mounting it onto a copper heating 
plate (see Figure 23). The copper heating plate would change the temperature of the 
sensor in a precise and incremental manner for the integrated thermal readout 
characterization. The temperature is swept from 297 K to 303 K (24° C to 30“ C) with 0.1 
K increments using a Thorlabs TC200 temperature controller. The light beam reflected 
off by the pixels over the CCD camera would deflect and be correspondingly fdtered by 
the 4F imaging system as the temperature increases. 

The change in intensity of the reflected light is monitored by the image processing 
software using background subtraction as described in earlier sections. An area on the 
sensor, corresponding to 3 x 3 pixels is selected to measure their average gray values as 
temperature changes. To ensure that the maximum thermal response is extracted from the 
test, the sensor FPA is also rotated on a rotating platform mount to compensate for the 
intrinsic bending of the sensor. The 3x3 pixels, which are selected to be tested, are based 
on initial test runs that showed they were totally dark at 24“ C and reflecting maximum 
amount of light at 30“ C. These pixels also need to be next to each other in the array. Each 
pixel’s gray level response is individually measured and the permanently dark areas of 
the substrate are not considered. The average values of the gray levels are then computed. 
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The average gray values from the 3x3 pixel area as a funetion of temperature ehange are 
shown in the graph in Figure 24. 
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Figure 24. Plot of average gray level as a function of temperature of the sensor 

array. 

From the data, it can be observed that the thermal response of the integrated 
system for the operating temperature range of 297 K to 303 K (24°C to 30‘’C) is 
approximately 23 gray levels per Kelvin. 
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c. 


EFFECT OF APERTURE SIZE 


The effect of the aperture size becomes apparent when imaging is performed 
using the integrated optical readout system. Using the same setup used for the system 
thermal response characterization, the background subtracted images for different 
aperture sizes are obtained for the same temperature range. 




(b) 

Figure 25. Background subtracted images for aperture size (a) 0.5 mm 

(b) 5 mm. 
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From Figure 26, it can be observed that aperture size dramatically changes the 
nature of the image produced by the optical readout system. It can be observed that the 
smaller aperture size results in the better filtering of lower spatial frequencies in the 
diffraction pattern in the Fourier plane. A larger aperture size results in the appearance of 
the higher spatial frequencies of the diffraction patterns. While the lower spatial 
frequencies are displaced as well, they are not filtered, and hence they do not show up in 
the background-subtracted image. This is evident in Figure 25(b) where the dark rings, 
which are the lower spatial frequencies, and the faint glow surround the dark rings are the 
higher spatial frequencies of the diffraction pattern. 


45 



THIS PAGE INTENTIONALLY LEET BLANK 


46 



VI. CONCLUSIONS 


In this thesis, an optical readout scheme was successfully developed based on the 
Fourier 4F optical system. The optical readout system was also successfully integrated 
with the MEMS THz detector. A QCL provided the THz illumination while post image 
processing performed background subtraction to obtain the THz scene on the detector. 

The concept of using optical readout to probe the detector was validated, as the 
integrated system was able to reproduce image of the QCL beam pattern due to the 
relative deflection of a large number of pixels. This was an important step in developing 
an image acquisition system based on optical readout. In addition, the physics of the 
Fourier 4F system was studied and image outputs due to spatial fdtering were checked 
and compared with established theories. This is to provide confidence in the ability of the 
Fourier 4F optical readout system to be able to enhance the sensitivity of the readout 
system. The sensitivity of the system is based on reflected light being blocked by the 
aperture as the sensor is deflected as well as optical spatial filtering. 

Through these studies, it was noted that there is room for the optical readout 
system to be further optimized for future work. For better collimation and lesser loss of 
THz radiation energy, a custom GAP mirror can be mounted on the QCL copper head to 
collimate and focus the THz laser before exiting the cryostat assembly. (This requirement 
was noted during the assembly of the optical system. However, the GAP was still in the 
process of fabrication and not available for testing.) This allows more energy to reach the 
MEMS THz detector to excite a greater number of pixels for better image reconstruction. 

In addition, the components used in the Lourier 4L optical system can be 
improved. The tested setup had the components aligned on an optical rail, which had 
considerable light leakage. Moving forward, lens tube systems can be designed around 
the Lourier 4L optical system concept to enhance the overall system. It is envisaged that 
lens tube systems can reduce distortions and aberrations as well as make the overall 
system more compact. 
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